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ABSTRACT 

The existence of a shallow decay phase in the early X-ray afterglows of gamma-ray bursts is a 
common feature. Here we investigate the possibility that this is connected to the formation of a 
highly magnetized millisecond pulsar, pumping energy into the fireball on timescales longer than 
the prompt emission. In this scenario the nascent neutron star could undergo a secular bar-mode 
instability, leading to gravitational wave losses which would affect the neutron star spin-down. In 
this case, nearby gamma-ray bursts with isotropic energies of the order of 10 50 ergs would produce a 
detectable gravitational wave signal emitted in association with an observed X-ray light-curve plateau, 
over relatively long timescales of minutes to about an hour. The peak amplitude of the gravitational 
wave signal would be delayed with respect to the gamma-ray burst trigger, offering gravitational 
wave interferometers such as the advanced LIGO and Virgo the challenging possibility of catching its 
signature on the fly. 

Subject headings: gamma rays: bursts; radiation mechanisms: non-thermal; gravitational waves 



1. INTRODUCTION 

Thanks to Swift obser vations (e.g. iNousek et al. I 
120061 IZhang et al.ll2006h . it has now become evident 
that the "normal" power-law behavior of long GRB X- 
ray light curves F(T) oc T a with a ~ —1.2 (where 
F(T) is the observed flux and T is the observer's time), 
is often preceded at early times by an initial steep de- 
cay (a ~ — 3), followed by a shallower-than-normal de- 
cay (a > —0.5, see Fig. Q}. The steep-to-shallow and 
shallow-to-normal decay transitions are separated by two 



corresponding break times, 100 s < Tf,. 



ik.l 



< 500 s 



and 10 3 s < Tb rea k,2 < 10 4 s. During the shallow-to- 
normal transition the spectral index does not change and 
the decay slope after the break (a ~ —1-2) is gener- 
ally consistent with the standard aftergl ow model (e.g. 
iMeszaros fc Rees Ill997t ISari et al. 1 119981 ). while the de- 
cay slope before the break is usually much shallower. The 
lack of spectral changes suggests that the shallow phase 
may be attributed to a continuous energy injection by 
a long-lived central engi ne, with progressiv ely reduced 
activity (for a review see IZhang et al. II2006I. and refer - 
ences therein). Recently. iPanaitescu fe Vestrand I (2008) 
have pointed out that the effects of a late-time energy 
injection may also be evident in some optical afterglows, 
around 30 — 10 4 s after the trigger. Although it is still 
not clear if a typical "steep-flat-steep" behavior does 
exist also in short GRB X-ray afterglows, the case of 
GRB 051221a does fit this scheme remarkably, with a 
plat eau observed right in th e middle of the afterglow de- 
cay (So derberg et. alll2006l ). 
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Newborn magnetars are among the progeni- 
tors proposed to account for shallow decays or 
plate a us observed in GRB ligh t curves dDai fe Lu I 
19M [Zhang fc Meszarosl 120011 : iFan fc Dong I 120061 : 
Yu fc Huang I 120071 ). Independent support for this 
scenario comes from the observation of SN2006aj, 
associated with the nearby sub-energetic GRB 060218, 
suggesting that the supernova-GRB connection may 
extend to a much broader range of stellar masses than 
previously thought, possibly involving two different 
mechanisms: a "collapsar" for the more massive stars 
collapsing to a black hole (BH), and a ne wborn neutron 
star (NS) for the less massive ones (M azzali et al. I 
120061) . 

Previous studies aimed at accounting for the after- 
glow plateaus by invoking a magnetar-like progenitor 
have assumed that the magnetar's slow-down is domi- 
nated by magnetic dipole losses, neglecting the contri- 
but ion from the emi s sion of gravitationa l waves (GWs, 
see iDai fe Lu I fl998t IFan fe Dong! 12001 lYu fe Huangl 
l20"07ir or treating them separately as a limiting case 
for a NS with sufficiently high, constant eccentricity 
([Zhang fe Meszaros I [2001 . Such studies have shown 
how magnetars dipole losses may indeed explain the flat- 
tening observed in GRB afterglows. In the simplest ver- 
sion of the magnetar scenario, the end of the shallow 
decay is accompanied by an achromatic break, while sev- 
eral cases of chromati c breaks have also been observed 
(e.g. IPanaitescu IfeOOQf ). Additional mechanisms such as 
variable m icro-physical parameters in the fireball shock 
front (e.g. IPanaitescu et al. 1 12006D or a structured jet 
model (e.g. iRacusin et al. ll200c^ U can be invoked to ex- 
plain such chromatic breaks. Anyhow, a larger sample 
of simultaneous optical-to- X-ray observations is needed 
to firmly asses the achromatic or chromatic behavior of 
breaks associated to the end of the shallow-decay phase. 

In this paper, we investigate in more detail the ef- 
fects of GW losses on the magnetar's spin-down, and 
explore quantitatively the signatures which could test 
whether this is indeed the mechanism at work in the 
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shallow X-ray light curves. Although the precise evolu- 
tion of a newborn magnetar from birth up to timescales 
of ~ 10 3 — 1 4 s is difficult to predict or to follow with 
numerical simulations, here we point out that among the 
possible evolutionary paths which one may reasonably 
consider, one plausible and particularly interesting pos- 
sibility to explore is that of a newborn NS left over after 
a GRB explosion, which undergoes a secular bar-mode 
instability. In this scenario, simple estimates accounting 
for the most relevant energy loss processes can provide 
useful insights into the viability of having efficient GW 
emission associated with a GRB X-ray afterglow plateau. 
Although these estimates are clearly approximate, since 
possible complications like viscosity effects or magnetic 
field driven instabilities are neglected, they nonetheless 
allow us to make a first statement on the relevance of the 
considered process. Moreover, while other scenarios are 
also possible, the interesting aspect of this particular one 
is that, on the one hand, GW observations would be facil- 
itated by the presence of an electromagnetic signature to 
pinpoint the GW signal search, while on the other hand 
the detection of bar- mode like GWs in coincidence with a 
GRB X-ray plateau would be a smoking-gun signature of 
a magnetar pumping energy into the fireball, thus identi- 
fying the much-debated plateau mechanism. Given that 
several alternative scenarios have been invoked to explain 
the afterglow flattening, which are no t expected to be as - 
sociated with GW signals (see e.g. iPanaitescu I [2008). 
this would represent a significant step forward in our un- 
derstanding of GRB physics. Moreover, identifying the 
presence of a magnetar would confirm that not all GRB 
explosions necessarily lead to the prompt formation of a 
BH. 

A point of interest for current analyses that GW de- 
tectors are carrying o ut (see e.g. lAcernese et al. I [2008; 
lAbbott et al. I l2008al lbl) is that the scenario described 
here involves a new class of GW signals, which should 
be searched for in coincidence with GRBs. These would 
have a longer duration (I0 3 — 10 4 s) and a different fre- 
quency evolution than the type of GW signals currently 
considered to be possibly associated with GRBs. More- 
over, being delayed by minutes to < 1 hour with respect 
to the prompt 7-ray trigger, the GW signal associated 
with a GRB plateau would offer the challenging possibil- 
ity of an on-line detection. In light of the fact that the 
Virgo 4 and LIGO 5 interferometers are now progressing 
toward their enhanced/advanced configurations, and get- 
ting prepared for performing on-line data analyses, this 
prospect appears very appealing. It is worth noting that 
despite a GW signal in coincidence with a GRB plateau 
could also be searched off-line by LIGO or Virgo, an on- 
line detection would be highly preferable, since it could 
serve as a trigger for ground-based optical follow-ups, 
even if a GRB trigger alert is absent for any reason. 

The paper is organized as follows. In Sec. [5] we briefly 
describe how GRB afterglow plateaus are modeled in the 
context of the magnetar model. In Sec. [3] we review 
the main processes that can lead to GW emission as- 
sociated with NS formation. The aim of this section is 
to show that, among the different mechanisms that can 
come into play, the high efficiency of the secular bar- 
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Fig. 1. — Cartoon representation of the typical light curve be- 
havior observed by Swift XRT. The "standard" power-law decay 
with index a = —1.2 is preceded by a flat phase, lasti ng 10 2 — 10 4 s, 
durin g which the decay index is a = —0.5 or flatter II Zhang et al. I 
l2006h . 



mode instability is conducive to producing GW signals 
which are detectable also from relatively nearby extra- 
galactic sources. Moreover, it develops on timescales 
compatible with the observed durations of GRB plateaus. 
Sec. [^describes the general idea and particular aspects of 
the scenario being explored here, and how it can explain 
GRB afterglow plateaus with the presence of a magne- 
tar whose spin-down includes both magnetic dipole and 
bar-mode GW losses. In Sec. [5] we present the results of 
our calculations, and in Sec. [H]we discuss these results, 
summarizing our conclusions in Sec. [7] 

2. GRB PLATEAUS IN THE MAGNETAR SCENARIO 

Although a wide range of GRB progenitors end 
in the formation of a BH-debris torus system, it 
has been proposed that some progenitors may lead 
to a highly magnetized rapidly ro t ating pulsar (e.; 
Usov I [19921: iDuncan fc Thompson I Il992t [T hompson 
199 



lYTfc Blackman I 119981: iBlackman fc Yi I 
Dai fc Lu~l 119981: iKluzniak fc Ruderman I 
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Nakamural Il998t [Spruit I Il999h IWheeler et all . 
Ruderman et al. 1 120001; iLevan et al. 1 120061 : iMereghetti I 
120081 : iBucciantini et al. I l2009h . with such possibility 
being realized not only in the case of long GRBs 
associated to collapsars, but eventually also in scenarios 
relevant for shor t GRBs, such as NS binary mergers 
([Dai fc Lu II1998L and refer ences therein). 

Fast rotating highly magnetized pulsars, are among 
the class of progenitors that may be associated with 
significant energy input in the fireball for timescales 
longer than the 7-ray emission, thus being relevant for 
explaining GRB afterglow plateaus. A detailed analy- 
sis of the observable effects linked to the presence of a 
pul sar pumping energy into t he fireball was performed 
by IZhang fc Meszaros I (|2001l ) , the results of which we 
briefly recall in what follows. 

Consider the general scenario where the GRB is pow- 
ered by a central engine that emits both an initial impul- 
sive energy input, Ei mp , as well as a continuous luminos- 
ity, the latter varying as a power-law with time, i.e. L — 

Lq ^Tjj) where T is the observer's time. This could be 

the case if the central engine is a pulsar and the initial im- 
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pulsive GRB fireball is due to v-v annihi lation or magne- 
tohydrodynamical processes (see e.g. iMacFadven et al~ 



f999; 



2005; 



Popham et al. Hl999tlDi Matteo et al. 
Oechslin fe Janka 1 120061 ; IZhang fe Dai 



2002tlLee 
2009]). In 



such self-similar blast wave is expected to form 

at late times. In a GRB, the timescale To at which 
the self-similar solution applies is roughly equal to the 
time for the external shock to start decelerating while 
collecting material from the interstellar medium (e.g. 
iSari fc Piran I [1999). At different times, the total en- 
ergy into the fireball may be dominated by either the ini- 
tial impulsive term, or by the continuous injection one, 

whose contribution will scale as Ei n j — L °^° {^f^j 

The continuous energy injection term can dominate on 
the impulsive one for T > T c (where T c > To so as to as- 
sure that the self-similar solution has already developed 
when the continuous injection law dominates), if q > —1 
and Ei n j(T c ) ~ Ei mp . In the particular case in which 
ToT) ~ Ei mp then T c ~ To and the dynamics is dom- 
inated by the continuous injection as soon as the self- 
similar evolution begins. Generally speaking, one can 

write T c = maxl T ,T [{q + l)E mp /(L T )] 1/{1+q} } 

(|Zhang fc Meszarosl |2001D . Note that the continuous 
injection may, in addition, have another characteristic 
timescale Tf at which the continuous injection power- 
law index q > — 1 switches to a lower value q < — 1. In 
such a case, it is only for T c < Tf that the continuous 
injecti on has a noticeable effect on the afterglow light 
curve ([Zhang fc Meszaros Il200lh . 

During the energy-injection dominated phase, the 
peak flux, peak frequency and cooling frequency of 
the sy nchrotron photon s produced by the forward 
shock ([Sari et al. I Il998| ) scale with time as F m cx 
T 1+ g, v m oc T-( 2 -g>/ 2 , y c oc T"(9+ 2 )/ 2 , respectively 
([Zhang fc Meszarosl 120011). that reduce to the standard 
scalings for q = — 1 ([Sari et al. lll998D . In the case of 
a nearly constant energy supply, i.e. q ~ 0, one has 
F m oc T, v m cx T _1 , v c cx T _1 , respectively. These 
scalings allow one to compute the temporal indices of 
the afterglow light curve expected during the injection 
phase. Supposing to be in slow cooling, these are F v cx 
T (3- P )/2 for Vm<v <Vc an d p v (x 

T^-P)/ 2 for v > v c , where we 
have indicated with p the power-law in dex of the electro n 
energy distribution in the shock front (|Sari et al. Ill998[ ). 
For 2 < p < 4, one has 0.5 > a% > —0.5 at frequencies 
v m < v < v c and > oli > — 1 at v > i> c , to compare 
with a > —0.5 observed during GRB afterglow plateaus. 
In the absence of energy injection, for the standard adi- 
abatic fireball one would have —3/4 > a± > —9/4 for 
v m < v < v c , and —1 > a>2 > —5/2 at v > v c , for the 
same range of p values. Thus, the presence of a pul- 
sar pumping energy into the fireball at a nearly constant 
rate, is expected to cause a flattening in the typical decay 
of the afterglow light curve, with a > —0.5, in agreement 
with Swift observations (see Fig. [1]). 

3. GWS BY NS FORMATION 

Gravitational collapse leading to the formation of a NS 
has long been considered an observable source of GWs. 
During the core collapse, an initial signal is expected to 
be emitted due to the changing axisymmetric quadrupole 



F r Jv- 1)l2 oc T ai 
F m 4r 1)7 V /2 oc T- 



moment. A second part of the GW signal is produced 
when gravitational collapse is halted by the stiffening 
of the equation of state above nuclear densities and the 
core bounces, driving an outwards moving shock, with 
the rapidly rotating proto-neutron star (PNS) oscillat- 
ing in its axisymmetric normal modes. In a rotating 
PNS, non-axisymmetric processes can also yield to the 
emission of GWs with high efficiency. Such processes are 
convection inside the PNS and in its surrounding hot en- 
velope, anisotropic neutrino emission, dynamical insta- 
bilities, and secular gravitational-radiation driven insta- 
bilities, that we br iefly recall i n wha t follows. We refer 
the reader to e.g. Kokkotas (2008) for a recent, more 
detailed review. 

- Convection a nd neutrino emission - 2D simulations 
of core collapse (Miil ler et al. 1 120041 ) have shown that 
the GW signal from convection significantly exceeds the 
core bounce signal for slowly rotating progenitors, being 
detectable with advanced LIGO for galactic sources. In 
many simulations, the GW signature of a nisotropic neu- 
trino emission has also been considered ([Epstein 1119781 ; 
iBurrows fc Haves Ill996t iMuller fc Jankalll997D and es- 
timated to be detectable by advanced LIGO for galactic 
sources. 

- Dynamical instabilities - They arise from non- 
axisymmetric perturbations and are of two different 
types: the classical bar-mode instability and the more 
recently discovered low-T/|W| bar-mode and one-armed 
spiral instabilities. In Newtonian stars, the classical 
m = 2 bar-mode instability is excited when the ratio 
(3 = T/\W\ of the rotational kinetic energy T to the grav- 
itational binding energ y \W\ is larger than j3d yn = 0.27 
( Chandrasekhar 1969). It can be excited in a hot PNS 
a few ms after core bounce, or alternatively a few tenths 
of seconds later, when the PNS cools due to neutrino 
emission and contracts further, with /3 becoming > (3d yn 
{j3 oc 1/R during contraction). The instability grows 
on a dynamical timescale (the time that a sound wave 
needs to travel across the star) which is about one ro- 
tational period, and may last from 1 to 100 rotations 
depending on the degree of differentia l rotation (e.g. 
iBaiotti et aT~ll2007t iManca et aT~ll2007ft . If the bar per- 
sists for ~ 10-100 rotation periods, then even signals from 
distances considerably larger than the Virgo Cluster are 
estimated to be detectable. An m = 1 one-armed spiral 
instability has also been shown to become unstable in 
PNS, provided that the differential rotation is sufficiently 
strong (with matter on the axis rotat ing at least ten 
times faster than matter on the equator. ICentrella et al. I 
120011 : [Saiio et al.~H 2002). In recent simulations of rotat- 
ing core collapse to which differential rotation was added 
(|Ott et al. II2005D , the emitted GW signal reached a max- 
imum amplitude comparable to the core-bounce axisym- 
metric signal, after ~ 100 ms and at a frequency of ~ 800 
Hz. 

- Secular instabilities - At lower rotation rates, a star 
can become unstable to secular non-axisymmetric in- 
stabilities, driven by gravitational radiation or viscos- 
ity. Secular GW-driven instabilities are frame-dragging 
instabilities usually called Chandrasekhar-Friedm an- 
Schutz (CFS. iChandrasekhar |[l970t iFriedman II1978D in- 
stabilities. Neglecting viscosity, the CFS-instability is 
generic in rotating stars for both polar and axial modes. 
In the Newtonian limit, the I — m — 2 /-mode, which 
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has the shortest growth time of all polar fluid modes 
(1 s < row < 7 x 10 4 s for 0.24 > p > 0.15, see 
lLai fc Shapiro Ill995t ). becomes unstable when (3 > 0.14. 
The /-mode instability, also referred to as the secu- 
lar bar-mode instability, is an ex cellent source of GWs . 
In the ellipsoidal approximation, lLai fc Shapiro I (|1995l ) 
have shown that the mode can grow to a large nonlin- 
ear amplitude, modifying the star from an axisymmet- 
ric shape to a rotating ellipsoid, that becomes a strong 
emitter of GWs until the star is slowed-down towards 
a stationary state. This stationary state is a Dedekind 
ellipsoid, i.e. a non-axisymmetric ellipsoid with inter- 
nal flows but with a stationary (non-radiating) shape 
in the inertial frame. During the evolution, the non- 
axisymmetric pattern radiates GWs sweeping through 
the advanced LIGO/Virgo sensitivity window (from 1 
kHz down to about 100 Hz), which could become de- 
tectable out to a distance of more t han 100 Mpc. Two 
recen t hydrodynamic al simulations ( Shi bata fc Karino I 
l2004UOu et al. |[2004l in the Newtonian limit and using 
a post-Newtonian radiation-reaction, respectively) have 
essentially confirmed this picture. 

Among axial modes, the I = m = 2 r-mode 

is an important member (see e.g. lAndersson I 

1998b iFriedman fc Morsink I 119981; iLindblom et all 
1998t lOwenetal. I 119981: iLindblom fc Owen I 120021 : 
Andersson fc Kokkotasl 120011: lAndersson I 120031: 
Bondarescu et al. I I2009D . If the compact object is 
a strange star, such instability is predicted to persist 
for a few hundred years (at a low amplitude) and, 
integrating data for a few weeks, could yield to an 
effective amplitude h c s ~ 10 ~ 21 for galactic signals, at 
frequencies - 700 - 1000 Hz (|Kokkotas 1120081) . 

- Other magnetic-field related effects - Finally, mech- 
anisms different from rotational instabilities can be in- 
voked as GW sources in newborn magnetars. E.g., 
in several scenarios the star's shape may be domi- 
nated by the distortion caused by very high internal 
magnet i c fiel d s (e.g. iPalomba] l2000l: I Cutler I 120021: 
Arons I l200l iStella e t al. I 12003: lDaH'Osso fc Stella I 
20071 iDairOsso et"al~l 12001) . GW signals produced by 
these kind of processes are typically estimated to be de- 
tectable by the advanced interferometers up to the Virgo 
Cluster (i.e. distances of the order of 20 Mpc). 

4. THE NS SPIN-DOWN 

On the longer afterglow timescales that are of interest 
for the present work, the energy injection into the fireball 
by a magnetar eventually surviving after the GRB explo- 
sion, is expected to be mainly through elec tromagnetic 
dipolar emission (|Zhang fe M eszaros! 1200 lh . For what 
concerns GW losses, in this work we focus on the secu- 
lar bar-mode instability, given its high efficiency in the 
production of GWs, and being its characteristic timescale 
tgw compatible with the one of GRB plateaus (see Sects. 
E]and[l]). 

As discussed in the previous section, a collapsing core 
rotating sufficiently fast is expected to become non- 
axisymmetric when (3 is sufficiently large. Since a new- 
born NS can be secularly unstable but dynamically stable 
only if the rotation rate of the pre-collapse core lies in a 
narrow range, and sin ce during the co l lapse (3 increases 
proportionally to i?~ 1 . lLai fc Shapiro! (|1995l ) considered 
more likely that the core becomes dynamically unstable 



(0 > Pdyn) following the collapse, provided the initial 
Pi is not too small. On a short dynamical timescale, 
such NS will evolve toward a nearly axisymmetric equi- 
librium state, with j3 decreas ing below /^n,, but p ossi- 
bly remaining above (3 sec (see lLai fe Shapiro 1119951 and 
references therein). Due to gravitational radiation, the 
nearly axisymmetric core (secularly unstable Maclaurin 
spheroid) will evolve into a non-axisymmetric configu- 
ration (Riemann-S ellipsoid), on a secular dissipation 
timescale ~ tqw- While an initial dynamical unsta- 
ble phase would possibly produce a GW burst during 
the GRB, the secular evolution takes place on longer 
timescales, thus being relevant for the shallow phase 
(100 s < T < 10 4 s) observed in GRB afterglows (see 
Fig. [J). For this reason, in what follows we focus on 
the secular bar-mode instability. It is worth noting, 
however, that also the presence of a bar-like GW burst 
from a dynamical bar-mode instability, would provide a 
hint for a magnetar being formed in the GRB explosion. 
BH formation, in fact, is not expected to lead to strong 
quadrupole moments (except i f it is argued for blobs 
formi ng in the infall, see e.g. iKobavashi fc Meszaros I 
2003), and in any case a dynamically unstable magne- 
tar would presumably give rise to a more regular signal. 

Fully general relativistic axisymmetric simulations of 
rotating stellar core collapse in three spatial dimension, 
performed for a wide variety of initial conditions (rota- 
tional velocity profile, equations of state, total mass), 
indicate that the threshold (3 = 0.27 for the onset of the 
classical dynamical instability is passed if the the progen- 
itor of the collapse is: (i) highly differentially rotating; 

(ii) moderately rapid ly rotating with 0.01 < Pi S, 0.02; 

(iii) massive enough (|Shibata fc Sekiguchi 1120051 ). More 
recent numerical collapse simulations of rotating stellar 
iron cores to PNS have also provided an extensive set of 
post-bounce rotational configurations, allowing studies of 
the prospects for the develo pment of non-axisymmetric 
rotational instabilities. E.g. iDimmelmeir et al~l (2008) 
found that the rotational barrier imposed by centrifugal 
forces prohibits the spin-up necessary for the classical dy- 
namical bar-mode instability, but a large subset of post- 
bounce models exhibits a p above t he secular instability 
thresh old. Based on their results, Dimmel meir et aLl 
(2008) consider it unlikely that a PNS in nature de- 
velops a high-/3 dynamical instability at or early after 
core bounce. While many of the PNS could theoreti- 
cally reach Pd yn during their subsequent cooling to the 
final condensed NS (if angular momentum is conserved), 
it is however considered more likely that the secular in- 
stability driven by dissipatio n or gravitational radiatio n 
back-reaction will set in first (jDimmelmeir et al. 1120081 ) . 
Still, three-dimensional simulations are necessary to pro- 
vide conclusive tests of these predictions. 

Under the hypothesis that a secular bar-mode insta- 
bility does indeed set in, in this work we follow the NS 
quasi-static evolution under the effect of gravitational ra- 
diation according to t he analytical formulation given by 
lLai fc Shapiro I (|1995f ). Such evolution can in principle 
be studied usi ng the full dynamical equations of ellip- 
soidal figures (|Chandrasekhar I IT969). including gravita- 
tional radiation reaction. However, since tqw is gener- 
ally much longer than the dynamical time of the star, 
the evolution is quasi-static, i.e. the star evolves along 
an equilibrium sequence of Ricmann-S ellipsoids. Dif- 
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ferently from what done bv lLai &; Shapiro! (|1995f) . here 
we add in the energy losses the contribution of magnetic 
dipole radiation, under the hypothesis that those will not 
substantially modify the dynamics, but will act speeding 
up the overall evolution of the bar along the same se- 
quence of Riemann-S ellipsoids that the NS would have 
followed in the absence of radiative losses. As we are 
going to show in the following section, dipole losses are 
nearly constant during the bar evolution so that, accord- 
ing to what discussed in Sec. [5] for the q = case, they 
can act as a source of continuous energy supply into the 
fireball, explaining the observed slope of GRB afterglow 
plateaus. 

It is worth noting that in a real situation where 
magnetic field instabilities and viscosity effects are also 
present, the relevant timescales may be altered. A secu- 
larly evolving bar can last up to a timescale of the order 
of 10 3 s, as far as viscosity or magnetic field induced in- 
stabilities do not substantially modify the dynamics. Vis- 
cosity may play a role on the secular evo lution when the 
PNS has cooled to below ~ 1 MeV (e.g. lLai fc Shapirol 
ll995tlLai"ll2001l ). The time required for the pulsar to cool 
below such te mper a ture w as estimated as a few hundreds 
of seconds by lLai I (|2001h . In the case of a GRB explo- 
sion, a less rapid cooling is expected due to continuous 
in-fall and j et emission (a heating source which was ab- 
sent in e.g. lLai Il200ll ). so that we may assume that the 
bar survives at least until the end of the electromagnetic 
plateau (i.e. T ~ 10 3 s). Magnetic effects are notori ously 
difficult to predict (see e.g. lShibata fc Karino Il200^ . and 
in general require making heuristic assu mptions. In the 
contex t of the secular r-mode instability, Rez zolla et al. I 
(2001) have shown that the growth of an initial magnetic 
field associated with the secular kinematic effects emerg- 
ing during the evolution of the instability, possibly damps 
the growth of the instability itself. Despite the different 
context (r- modes) , these results do suggest that magneti- 
cally driven instabilities may complicate the scenario. In 
what follows, we explore the quantitative consequences of 
making the plausible assumption that magnetic instabil- 
ities are less efficient at spinning-down the bar than GW 
emission and magnetic dipole losses. Moreover, apart 
from magnetic braking (spin-down due to dipole losses) 
which we do consider here, the presence of a magnetic 
field can influence the secular evolution in other ways. A 
magnetic field anchored on the star's surface is in fact 
perturbed by the instability itself, and this can lead to 
electromagnetic losses which can enha nce the CFS mech- 
anism. E.g. in the context of r-modes. lHo fc Lai I (2000) 
have considered the electromagnetic radiation associated 
with the shaking of magnetic field lines by the r-mode 
oscillations. This effect has been estimated to be negli- 
gible for NS with magnetic field strengths below 10 16 G. 
In our case, magnetic field lines anchored on the surface 
would be distorted by the bar-mode instability. In our 
treatment, we neglect this effect, and include only dipole 
losses associated with a magnetic field whose flux is con- 
served on a sphere of radius equal to the mean radius of 
the ellipsoid. 

5. MODELING OF THE NS EVOLUTION 

A general Riemann-S ellipsoid is characterized by an 
angular velocity fi e~3 of the ellipsoidal figure (the pattern 
speed) about a principal axis £3, and by internal fluid 



motions which are assumed to have uniform vorticity £ e$ 
along the same axis (in the frame co-rotating with the 
figure). Labeling with a\ and a 2 the principal axes of 
the ellipsoidal figure in the equatorial plane 6 , and with 
x\ and x 2 Cartesian coordinates in such a plane, it can be 
shown that the fluid velocity in the inertial frame reads 



with 



uq = u + fi(e3 x r) 



u = — Ax2ei Ax\e2 

a 2 at 



1 



(2) 



the velocity in the frame co-rotating with the figure 
( Chandrasekhar 1969), where ei and e 2 are unit vectors 
along the Cartesian axes x\ and x 2) r is the position 
vector; x indicates the vector product; and 



J1 -2 , _2^ 



(3) 



is the angular frequency of the internal fluid motions, 
i.e. of the elliptical orbits that the particles span around 
the rotational axis in addition to the pattern motion. 
The velocity uq is contained in the plain perpendic- 
ular to the rotational axis e^, so indicating with r*j_ 
the projection of the position vector r in such plane 
we can write r — rj_ + £363, and uq = d(f±)/dt + 
{dx$/d£)e§ — d(f±)/dt (i.e. dx^/dt — since the compo- 
nent of uq along S3 is null). Further, we can write uq — 

(dr±/dt)rj_ + rj_d(r±)/dt = (dr±/dt)r± + r±(Q x fx), 
where we have set 



fin 



1 



(Pj_ X U ). 



(4) 



Note that 51 is defined in such a way that r±fi gives 
the component of the particle velocity perpendicular to 
the polar radius fj_, measured in the inertial frame. How- 
ever, as underlined above, the motion of fluid particles 
on the surface can be viewed as the superposition of a 
circular motion with the pattern frequency fi, plus an 
elliptical motion on paths contained on the pattern ellip- 
soid (resulting in maintaining the pattern fixed). Since 
the internal fluid motions are ellipses rather than circles, 
there is an additional component of the velocity parallel 
to r±. Using Eq. JT]), one has 



fin 



fi 



A 



e.3 



(5) 



In the frame co-rotating with the pattern, fluid par- 
ticles on the star's surface move around the rotational 
axis, on ellipses contained in x 3 = const planes. Those 
ellipses are self-similar to the equatorial one and have 
equation: 



(6) 



6 In the ellipsoidal approximation, surfaces of constant density 
are assumed to be self-similar ellipsoids, so the geometry of the 
configuration is completely specified by the three principal axes 
oi, a 2 and (13, and the axis rati o as/0.1 and (12/^1 are the sa me for 
all interior isodensity surfaces (Lai, Rasio & Shapiro 1993). 
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For such fluid particles, a\x\ + a\x\ = afa 2 , ( 1 — ■ 
and (r_i_) = Wflia 2 ( 1 - : 



so that in one cycle 



a 2 x 2 +a 2 x 2 



-. , — 1. Thus, in the inertial frame, fluid parti- 

r 2 L a 1 a 2 / ' ' 1 

cles on the star's surface are characterized by an angular 
frequency: 



a eff = (n ) = q-a 



(7) 



Since the gravitational radiation reaction acts like a 
potential force, the fluid circulation along the equator of 
the star, 



C = 



u - dl — naia 2 (o, 



(8) 



equator 



where dl is taken along the star's equator and Co is the 
vorticity in the inertial frame, is conserve d in the absence 
of viscosity ([Lai. Rasio fe Shapiro lfl993T ). Therefore, the 
NS will follow a sequence of Riemann-S ellipsoids with 
constant circulation. Treating the NS as a polytrope of 
index n (|Chand rasckhar 1939), total mass M, and in- 
dicating with R the radius of the non-rotating, spher- 
ical equilibrium polyt rope with same mass M, one has 
(|Lai fc Shapiro II1995D : 



C 



Mk n C 



(9) 



where G is the gravitational constant; k n is a con- 
stant which depend s on the index n of the con sidered 
polytrope (see e.g. lLai. Rasio fe Shapiro 1119 9 3h - Note 
that C = C/VGA/ 3 i?o is an adimensional quantity, C — 
— (k n MC) / (hit) has the dimensions of an angular mo- 
mentum, and both are proportion al to the conserved cir - 
culation C. It can be shown that (|Lai fc Shapiro Ill995l) : 



C = I A k n Maia 2 Sl 

5 



(10) 



where / = k n M(a\ -{-a 2 ,)/ 5 is the NS moment of inertia 
with respect to the rotational axis. Along the secular 
equilibrium sequence, we wr ite the NS spin-down law as 
([Shapiro fc Teukolskv lfl983h : 



dE 
dT 



6c 3 



32GI 2 e 2 il 6 
5c5 



L 



dip 



(ii) 



where E is the NS total energy, Law — dEcw/dT ac- 
counts for GW losses, while Ldi P = dEd ip /dT for mag- 
netic dipole ones. Here e = (a\ — al)/(a 2 + a l) is 
the ellipticity; B p is the dipolar field strength at the 
poles; n is the pattern angular frequency of the ellip- 
soidal figure; R is the mean stellar radius; c is the 
light speed; T is the time measured in an inertial frame 
where the pulsar is at rest. Ldip is computed conserving 
the magnetic field flux over a sphere of radius equal to 
the mean stellar radius (i.e. B p R 2 — const — B p flR^ 
along the sequence, where R is the geometrical mean of 
the ellipsoid principal axes), and using the effective an- 
gular frequency ^ e //, which includes both the pattern 
speed and the effects of the internal fluid motions. The 



use of Vl 



(Ho) 



\r± x 



accounts for the 



fact that in the frozen-in magnetic field approximation 
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□ 
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Fig. 2.— NS evolution along a Riemann-S sequence with pa- 
rameters (C,n,M,Ro,B P}0 ) = (-0.41,1,1.4 M Q , 20 km, 10 14 G). 
Upper panel: Rate of energy loss in units of 10 45 ergs/s, when both 
magnetic dipole losses (red-dashed line) and GW losses (black- 
solid line) are taken into account in the magnetar's spin-down 
law. For reference, we also plot the rate of energy loss in the 
ca se only GW em i ssion is considered (black-dash-dotted line), as 
in lLai fc Shapiro (1995|). Lower panel: absolute value of the sur- 
face fluid particles angular frequency divided by a factor of n (i.e. 
\Q e f f\/ir), when both magnetic dipole and GW losses are consid- 
ered (black-solid line). For reference, we also plot the same quan- 
tity when only GW losses are taken into account in the magne- 
tar's s pin-down law (black-dash-dotted line), as in Lai & Shapiro 
(1995). Note that the vertical axis in the lower panel is a linear 
scale: between 10 2 s and ~ 10 3 s, fl e ff/TT changes from ~ 800 Hz 
to ~ 750 Hz, i.e. less than ~ 10% of its initial value. Thus, be- 
tween 10 2 s and 10 3 s the power-law approximation to dipole losses 



is L dip oc T 



so that q ~ can be assumed for T < 10 s. 



(See the electronic version for colours). 



see e.g. iGoldreich fc Julian" 119691: iBavm et al. 



1969 



Thompson fc Duncan I [19961 : iMorsink fc Rezania 112002 : 
Thompson et al. 1 l2002f) the magnetic field lines are in 
effect tied to the fluid particles on the stellar surface. 
Note that fi e / / (and the corresponding dipole loss term) 
is measured in the inertial frame, which is where we com- 
pute dE/dT as well. 

Once (C,n, M, Rq, Bp.o) are assigned, each configura- 
tion along a constant-C sequence is completely deter- 
mined specifying the axis-ratio x = a 2 /ai in the ellip- 
soid equatorial plane. Thus, all relevant quantities can 
be considered as functions of x only, and Eq. pip can 
be written as: 



dx 
dT 



Ldi P {x) 



->GW 



(x) 



dE/dx 



(12) 



We solve the above equation numerically, with its right 
hand side evaluated along a constant-C Riemann-S se- 
quence, and imposing an initial condition sufficiently 
near to a uniformly rotating Maclaurin spheroid, (x(ti) — 
Xi — > 1) of the given circulation C. 

6. RESULTS AND DISCUSSION 

In Fig. [2] we compare the luminosity emitted in GWs, 
computed with (black-solid line) or without (black-dash- 
dotted line) the addition of the dipole loss term (red- 
dashed line) in Eq. (fTTj) . for a typical choice of pa- 
rameters, (C, n, M, i? _, B p<0 ) = (-0.41, 1, 1.4 M Q , 20 km, 
10 14 G). Note that C = —0.41 corresponds to a value 
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of (3 = 0.20 for the initial Maclaurin configuration, i.e. 
in the middle of the 0.14 < < 0.27 range for the 
secular instability. As evident from the lower panel of 
Fig. [2 as long as the circulation is conserved, fl e ff 
remains nearly constant during all the evolution, and 
\Ldi P \ ~ 3 x 10 47 ergs/s = L (upper panel, red-dashed 
line). As underlined in Sec. [2J energy pumped into the 
fireball at a constant rate is sufficient to explain the 
observed temporal behavior of afterglow plateaus (i.e. 
a > —0.5, see Fig. [T]). For what concerns the duration 
of the plateau, for a GRB with impulsive isotropic energy 
of the order of E imp ~ 10 50 ergs, the effect of the energy 
injection in the light curve will become visible after a 
time T c ~ E imp /L ~ (10 50 ergs)/(3 x 10 47 ergs/s) - 300 
s (see the red-dashed line in Fig. [2] and Sec. [3]), 
which is about in the middle of the observed range for 
Tbreak,i ~ 100 — 500 s (see Fig. [T]). Supposing the energy 
injection ends or starts fading significantly when the star 
approaches the final Dedekind state (see the discussion 
at the end of this section), the GRB light curve will re- 
turn to its standard behavior after Tb re ak,2 ^ 10 3 , to be 
compared with the observed range of 10 3 — 10 4 s. Thus, 
for a GRB with such impulsive energy, the properties of 
the plateau associated to the NS secular evolution are in 
agreement with those typically observed . 

The waveform of the GW signal emitted in as- 
sociation with the a fterglow plateau is computed as 
(|Lai fc Shapiro I [l995h : 



hit) , , 
h + = Y~ cos $(1 + cos i 



h x = —h(t) sin $ cos 9 
(13) 

where 9 is the angle between the line of sight and the 
rotation axis of the star, $ = 2 f. fit is twice the orbital 
phase, and 



h(t) 



'2c 3 d 2 & 
5G 



->GW\ 



c 4 d 



Ie 



(14) 



where d is the distance to the source, Lqw and f2 are 
shown in Fig. [2] (upper panel, black-solid line) and Fig. [3] 
(lower panel, black solid line), respectively. The resulting 
GW signal is quasi-periodic, with frequency / = Q,/it. 

To estimate the GW signal detectability, we proceed 
as follows. For broad-band interferometers such as LIGO 
and VIRGO, the best signal-to- noise ratio is obtained by 
applying a matched filtering technique to the data, when 
a waveform template is available. In such a case, 



P' 



= 4 



o 



{Fl\h + {fM 2 + Fl\hAJM 2 ) 
S h (f) 



df (15) 



where h is a Fourier transform; Sh(f) is the power 
spectral density of the detector noise; F + , F x are 
the beam pattern functions (0 < F 2 ,F 2 < 1 de- 
pending on the source position in the sky, see e.g. 
lThorne l[l987l iFlanagan fe Hughes I [1991 . For the sig- 
nal in Eg. ()13j) , in the stationary phase approximation 
(iThqrne I Il987t ICutler fe Flanagan I 11994 lOwen et al. I 
l!998H0wen fe Lindblom 112002ft . 

7 Note that in our discussion we are neglecting redshift effects, 
since we are interested in nearby GRBs at z < 0.035, i.e. having 
luminosity distances d^ < 150 Mpc. 
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Fig. 3. — Upper panel: characteristic GW amplitude h c at 
d = 100 Mpc, with dipole plus G W (black-solid line) a nd only GW 
(black-dash-dotted line, see also [Lai & Shapiro 1995) losses being 
considered. A typical fit to the sensiti vity expected for advanced 
detectors (purple-d ashed line, see e.g. ICut ler & Flanagan] 11994 : 
lOwen et al. I [19981) . Virgo nominal sensitivity (blue-dotted line), 
and the advanced Virgo sensitivi ty optimized for binar y searches 
(blue-dash-dot-dot-dotted line, Accrnesc et al. 2008), are also 
shown. Lower panel: evolution of the GW signal frequency, with 
dipole plus GW (black-solid line) and only GW (black-dash-dotted 
line) losses being considered in the NS spin down. (See the elec- 
tronic version for colours). 



„ f + °° h 2 (t) (dt/ df ■) (Fi (1 + cos 2 9) 2 /4 + F$ cos 2 9) 

p =h m df - 



(16) 



Since we expect to be observing the GRB on-axis, 9 
0. In case of optimal orientation, 



^ = / + "/™*), (ln/) = 

Jo Jbh{j) 



d(ln/), (17) 



being h c — fh(t)wdt/df the characteristic ampli- 
tude, and h rms = y/fSfJj). In the upper panel of 
Fig. [31 we compare h c computed for a GRB at d = 
100 Mp c, with the h rm s expected for the advanced de- 
tectors (lAcemese et al~ll2008l : ICutler fc Flanagan Ill994h 
lOwen et al. I [1998), for which p ma x ^ 5 at d < 
100 Mpc, or d < 150 Mpc if we make the assump- 
tion that knowledge of the GRB trigger time reduces 
the detection threshold, of a factor which as a rule-of - 
thumb we take e qual t o 1.5 (jKochanek fc Piranl 119931 : 
ICutler fc Thorne 1 12002ft . Higher confidence m an even- 
tual detection may require {p) sky = ^/2/5 p ma x ^ 5 
in each of a three-detecto rs network with similar h rms 
(|Cutler fc Flanagan I [l994ft . With the help of the GRB 
trigger time roughly compensating the factor of 
this implies d < 100 Mpc. 

BATSE results show that about 3% of short GRB s 
are expected to be within 100 Mpc (Na kar et al. 1 12006') . 
which translates into ^1 — 2 short GRBs per year 
in the Swift (~ 10 short GRBs per year) plus GBM 
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(GLAST - - 1/4 of - 200 GRBs per year) sample. As 
far as low- luminosity long GRBs, two of them (980425 
and 060218) were already observed at a d ~ 40 Mpc 
and d ~ 130 Mpc, and their local rate (> 200 Gpc -3 
yr" 1 ) is expected to be much hig her than that of nor - 
mal bursts (1 Gpc" 3 yr" 1 , e.g. IVirgili et al~l 12009( 1. 
INTEGRAL has detected a large proportio n of faint 
GRBs inferred to be local ((Foley et al. ll2008T l, a sample 
which may be increased b y future missions such as Janus 
(|Stamatikos et al. |[2009h and EXIST 8 . To compare with 
the case discussed here, the standard progenitor scenario 
for long GRBs predicts p ~ 5 at 27 Mpc for the advanced 
Virgo/LIGO, while the chirp signal from short GRBs is 
esti mated to be detectable up to several hundreds Mpc 
(e.g. iFlanagan fc Hughes 1998; Kob avashi fc Meszaros I 
2003). GWs eventually detected after a chirp and dur- 
ing an electromagnetic plateau of a short GRB, would 
add a significant piece of information, probing whether 
a magnetar is formed in the coalescence, rather than a 
BH. 

It is finally worth adding some few considerations 
more. First, in the scenario we are proposing here, some 
correlations do exist between the electromagnetic plateau 
and the GW signal, that could be explored in future anal- 
yses so as to test up to which level those may help the 
GW signal search. For example, a measurement of the 
initial frequency of the GW signal, for a given NS mass 
and radius, would allow one to derive a n estimate 0qw 
for th e actual value of (3 (see also Fig. 5 in lLai fc Shapiro I 
11995( 1 . In the ellipsoidal approximation, (3qw would pre- 
dict a specific evolution of the bar, as e.g. the expected 
value of fl e ff((3 GW ) = H(/3 G w) ~ HPgw) during the 
nearly constant phase. At the same time, the luminos- 
ity of the afterglow plateau, for a given NS mass, ra- 
dius and magnetic field strength, would also allow one 
to estimate the value of £l e ff during the constant phase, 
which could thus be checked for consistency with the 
value Q, e f f(pcw) inferred from the GW measurements. 

Next, some considerations are required on the fate of 
the bar after the final Dedekind state is reached. In the 
absence of dipole losses, the evolution of the NS along the 
sequence would have maintained f2 nearly constant up to 
a time Tqw, of the order of few secular growth times, 
tow - 2 x IP" 5 sJM/[1AM q)]~ 3 [i?o/(10km)] 4 [fi - 
[3 sec ) 5 (|Lai fc Shapiro I (199511 . Here [3 is referred to the 
initial Maclaurin configuration^ and it's determined by 
the choice of C. In our case, C = —0.41 and — 0.20, 
so that tqw — 335 s. As evident from the black-dash- 
dotted line in the lower panel of Fig. (3( for such value of 
the circulation, when only GW losses are considered, one 
has Tqw ~ (1 - 2) x 10 3 s ~ (3 - 6) r. The addition of 
magnetic dipole losses speeds-up the process, so that the 
star reaches the stationary football configuration some- 
what earlier (Fig. 02 lower panel, black-solid line). After 
the end of the secular e volution, we do not k now what 
the fate of the bar is. As lLai &: Shapiro! (|1995f ) have un- 
derlined, while the star approaches a Dedekind ellipsoid 
the gravitational evolution timescale increases, eventu- 
ally becoming comparable to the viscous dissipation one. 
When this happens, C is not conserved anymore and the 
star is expected to be driven along a nearly-Dedekind 
sequence to become a Maclaurin spheroid, since this is 
the only final state that does not radiate GWs or dissi- 



pate energy viscously. The addition of magnetic dipole 
losses would speed up such evolution, and further spin- 
down the final Maclaurin state. We thus expect to_have 
ri e // decreasing at some point after the constant— C evo- 
lution, with the dipole luminosity Ldip also decreasing 
accordingly. Correspondingly, the energy injected into 
the fireball will start decreasing (eventually entering in 
the q < — 1 phase, see Sec. [2]), and the afterglow plateau 
is expected to end, with the light curve turning back to 
the temporal decay expected in the absence of contin- 
uous energy injection. In view of these considerations, 
and for the purpose of this paper, we have limited our 
discussion to show that the properties of the electromag- 
netic plateau, even assuming this suddenly ends after the 
constant— C evolution, are in agreement with those typi- 
cally observed in GRBs. 

7. CONCLUSION 

We have discussed a possible scenario where a newly 
formed magnetar is left over after a GRB explosion, and 
explored the hypothesis of its being subject to a secu- 
lar bar-mode instability, including in the spin down the 
contributions of both radiative losses by magnetic dipole 
emissio n and by GWs. Follow ing the analytical treat- 
ment of lLai fc Shapiro I ((1995(1 . we have shown that for 
reasonable values of the physical parameters, the typical 
properties of GRB afterglow X-ray plateaus may be re- 
produced. A consequence of this is that, on the relatively 
long timescale 10 3 — 10 4 s of the electromagnetic plateau, 
the advanced LIGO/ Virgo interferometers may detect a 
corresponding GW signal up to d ~ 100 Mpc, by carry- 
ing out matched searches. Such a signal would be asso- 
ciated with an afterglow light curve plateau from a long 
sub-luminous GRB, or from a short GRB, with isotropic 
energy < 10 50 erg, which is typical of most nearby GRBs 
detected. For the more energetic GRBs, a bar-mode GW 
signal may be detected without a visible plateau in the 
afterglow. 

In conclusion, although there are considerable uncer- 
tainties about the evolutionary path of newborn mag- 
netars, our analysis indicates that the scenario pro- 
posed here is a plausible and interesting possibility, lead- 
ing to an efficient GW emission process which is ac- 
companied by a distinctive electromagnetic signature. 
Thus, in view of the impending commissioning of the 
advanced LIGO and Virgo, we consider that it would 
be highly worthwhile to test this possibility through 
matched electromagnetic- GW data searches. 



We are grateful to Benjamin Owen for important com- 
ments and valuable suggestions on this scenario, and for 
helping improve the manuscript. AC thanks Fulvio Ricci 
for crucial support during this project, Giovanni Mon- 
tani for important discussions, Cristiano Palomba for 
very helpful suggestions, and Christian D. Ott for use- 
ful comments. This work was supported by the "Fon- 
dazione Angelo della Riccia" - bando A. A. 2007-2008 & 
2008-2009 (AC), and by NSF PHY-0757155 & NASA 
NNX08AL40G (PM). AC gratefully acknowledges the 
support of the Penn State Institute for Gravitation and 
the Cosmos (IGC). 



GRB afterglow plateaus and Gravitational Waves 



9 



http:/ /exist. gsfc.nasa.gov/ 

REFERENCES 



Abbott, B. et al. 2008, Phys. Rev. D, 77, 062004 
Abbott, B. et al. 2008, ApJ, 681, 1419 

Acernese, F. et al. 2008, Class. Quantum Grav., 25, 225001 
Acernese, F. et al. 2008, "AdV Prelimiary Design", VIR-089A-08 
Andersson, N. 1998, ApJ, 502, 708 

Andersson, N., and Kokkotas, K. D. 2001, Int. J. Modern Phys. D, 
10, 381 

Andersson, N. 2003, Class. Quantum Grav., 20, 105 
Arons, J. 2003, ApJ, 589,871 

Baiotti, L., et al. 2007, Phys. Rev. D, 75, 044023 
Baym, G., et al. 1969, Nature, 224, 674 
Bondarescu, R. et al. 2009, Phys. Rev. D 79, 104003 
Blackman, E. G., and Yi, I. 1998, ApJ, 498, L31 
Bucciantini, N. et al. 2009, MNRAS, 396, 2038 
Burrows, A., and Hayes, J. 1996, Phys. Rev. Lett., 76, 352 
Centrella, J. M. et al. 2001 ApJ, 550, L193 

Chandrasekhar, S. 1939, "An Introduction to the study of Stellar 

Structure", Chicago: Chicago Univ. Press) 
Chandrasekhar, S. 1969, "Ellipsoidal Figures of Equilibrium" (New 

Haven: Yale Univ. Press) 
Chandrasekhar, S. 1970, Phys. Rev. Lett., 24, 611 
Culter, C, & Flanagan, E. E. 1994, Phys. Rev. D, 49, 2658 
Cutler, C. 2002, Phys. Rev. D 66, 084025 

Cutler, C, & Thorne, K. S. 2002, 16th Int. Conf. on Gen. Rel. and 
Grav., Ed: Bish op, N. T. and Sunil, D. M. (Singapore: World 
Scientific), p.72, |arXiv:gr-qc/0204090| 
Dai, Z. G., & Lu, T. 1998, A&A, 333, 87 
DalFOsso, S., & Stella, L. 2007, Ap&SS, 308, 119 
DalFOsso, S., Shore, S. N., & Stella, L. 2008, to appear in MNRAS 

(arXiv0811.4311D) 
Di Matteo, T., Perna, R., & Narayan, R. 2002, ApJ, 579, 706 
Dimmelmeir, H. 2008, Phys. Rev. D, 78, 064056 
Duncan, R. C, & Thompson, C. 1992, ApJ, 392, L9 
Epstein, R. 1978, ApJ, 223, 1037 
Fan, Y.-Z., & Dong, X. 2006, MNRAS, 372, 19 
Flanagan, E. E., & Hughes, S. A. 1998, Phys. Rev. D, 57, 4535 
Foley, S., et al. 2008, A&A, 484, 143 
Friedman, J. L., & Schutz, B. F. 1978, ApJ, 222, 281 
Friedman, J. L., & Morsink, S. 1998, ApJ, 502, 714 
Goldreich, P., & Julian, W. H. 1969, ApJ, 157, 869. 
Ho, W. C. H., & Lai, D. 2000, ApJ, 543, 386 
Kluzniak, W., & Ruderman, M. 1998, ApJ, 505, L113 
Kobayashi, S., & Meszaros, P. 2003, ApJ, 589, 861 
Kochanek, C. S., & Piran, T. 1993, ApJ, 417, L17 
Kokkotas, K. D. 2008, Rev. in Modern Astrophys., 20, 140 
Lai, D., Rasio, F. A., & Shapiro, S. L. 1993, ApJS, 88, 205 
Lai, D., & Shapiro, S. L. 1995, ApJ, 442, 259 

Lai, D. 2001, in AIP Conf. Proc. 575, Astrophysical Sources for 
Ground-based Gravitational Wave Detectors, ed. J. M. Centrella 
(Melville, NY: AIP), 246 
Lee, W. H., et al. 2005, ApJ, 632, 421 
Levan, A. J., et al. 2006, MNRAS, 368, 1 
Lindblom, L., et al. 1998, Phys. Rev. Lett., 80, 4843 
Lindblom, L., & Owen, B. J. 2002, Phys. Rev. D, 65, 063006 
MacFadyen, A. I., & Woosley, S. E. 1999, ApJ, 524, 262 



Manca, G. M., et al. 2007, Class. Quantum Grav., 24, S171 

Mazzali, P. 2006, Nature, 442, 1018 

Mereghetti, S. 2008, A&A Rev., 15, 225 

Meszaros, P., & Rees, M. J. 1997, ApJ, 476, 232 

Morsink, S. M., & Rezania, V. 2002, ApJ, 574, 908 

Miiller, E., & Janka, H. T. 1997 A&A, 317, 140 

Miiller, E., et al. 2004, ApJ, 603, 221 

Nakamura, T. 1998, Prog. Theor. Phys., 100, 921 

Nakar, E., et al. 2006, ApJ, 650, 281 

Nousek, J. A., et al. 2006, ApJ, 642, 389 

Oechslin, R., & Janka, H.-T. 2006, MNRAS, 368, 1489 

Ott, C. D. et al. 2005, ApJ, 625, L119 

Ou S. et al. 2004, ApJ, 617, 490 

Owen, B. J. et al. 1998, Phys. Rev. D, 58, 4020 

Owen, B. J. & Lindblom, L. 2002, Class. Quantum Grav., 19, 12470 

Palomba, C. 2000, A&A, 354, 163 

Panaitescu, A. et al. 2006, MNRAS, 369, 2059 

Panaitescu, A. 2008, MNRAS, 383, 1143 

Panaitescu, A., & Vestrand, W. T. 2008, MNRAS, 387, 497 

Panaitescu, A. 2009, in AIP Conf. Ser. 1133, Gamma-ray Burst: 

Sixth Huntsville Symp., ed. C. Meegan, C. Kouveliotou, & N. 

Gehrels (Melville, NY: AIP), 127 
Popham, R. et al. 1999, ApJ, 518, 356 
Racusin, J. et al. 2008, Nature, 455, 183 
Rezzolla, L., et al. 2001, Phys. Rev. D, 64, 104014 
Ruderman, M., Tao, L., & Kluzniak, W. 2000, ApJ, 542, 243 
Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 
Sari, R., & Piran, T. 1999, ApJ, 520, 641 
Saijo M. et al. 2002, ApJ, 595, 352 

Shapiro, S. L., & Teukolsky, S. A. 1983, "Black Hole, White Dwarfs, 
and Neutron Stars: The Physics of Compact Objects" (New 
York: Wiley) 

Shibata, M., & Karino, S. 2004, Phys. Rev. D, 70, 4022S 
Shibata, M. & Sekiguchi, Y. 2005, Phys. Rev. D, 71, 024014 
Soderberg, A. M. et al. 2006, ApJ, 650, 261 
Spruit, H. C. 1999, A&A, 341, LI 

Stamatikos, M. et al. 2009, Contributed to the SSE-SFP of the 

Astro2010 Decadal Survey (arXiv0902.3022S) 
Stella, L. et al. 2005, ApJ, 634, L165 
Thompson, C. 1994, MNRAS, 270, 480 
Thompson, C, & Duncan, R. C. 1996, ApJ, 473, 322 
Thompson, C, et al. 2002, 574, 332 

Thorne, K. S. 1987, "300 Years of Gravitation" , p. 330 (Cambridge: 

Cambridge Univ. Press) 
Usov, V. V. 1992, Nature, 357, 472 
Virgili, F. J. et al. 2009, MNRAS, 392, 91 
Wheeler, J. C, et al. 2000, ApJ, 537, 810 
Yi, I., & Blackman, E. G. 1998, ApJ, 494, L163 
Yu, Y., & Huang, Y.-F. 2007, ChJAA, 7, 669 
Zhang, B., & Meszaros, P. 2001, ApJ, 552, 35 
Zhang, B. et al. 2006, ApJ, 642, 354 

Zhang, D., & Dai Z. G. 2009, ApJ, submitted <|arXiv:0901.043H /l') 



